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Measuring the transport of electrons through a graphene sheet necessarily involves contacting 
it with metal electrodes. We study the adsorption of graphene on metal substrates using first- 
principles calculations at the level of density functional theory. The bonding of graphene to Al, 
Ag, Cu, Au and Pt(lll) surfaces is so weak that its unique "ultrarelativistic" electronic structure 
is preserved. The interaction does, however, lead to a charge transfer that shifts the Fermi level 
by up to 0.5 eV with respect to the conical points. The crossover from p-type to n-type doping 
occurs for a metal with a work function ~ 5.4 eV, a value much larger than the work function of 
free-standing graphene, 4.5 eV. We develop a simple analytical model that describes the Fermi level 
shift in graphene in terms of the metal substrate work function. Graphene interacts with and binds 
more strongly to Co, Ni, Pd and Ti. This chemisorption involves hybridization between graphene 
Pz-states and metal d-states that opens a band gap in graphene. The graphene work function is 
as a result reduced considerably. In a current-in-plane device geometry this should lead to n-type 
doping of graphene. 

PACS numbers: 73.63.-b, 73.20.Hb, 73.40.Ns, 81.05.Uw 



I. INTRODUCTION 

The history of carbon-based electronics begins with 
the discovery of fullerenes and carbon nanotubes that 
are zero- and one-dimensional, respective ly.^'^! A very 
important recent development was the preparation of 
single monolayers of graphite, now more commonly 
called graphene, on insulating substrates using mi- 
cromechanical cleavag^'^ that has made possible electron 
transp ort exper iments on this purely two-dimensional 
system.HISnillMl These transport measurements reveal 
high charge carrier mobilities, quantization of the con- 
ductivity, and a zero-energy anomaly in the quan tum 
Hall effect, as predicted theoretically.ESEHilEMl The 
theoretical studies explain these spectacular effects in 
terms of graphene's unique electronic structure. Al- 
though a single graphene sheet is a zero-gap semiconduc- 
tor with a vanishing density of states at the Fermi energy, 
it shows metallic behavior due to topological singulari- 
ties at the ii'-points in the Brillouin zone,^" where the 
conduction and valence bands touch in so-called coni- 
cal or Dirac points and the dispersion is essentially lin- 
ear within ±1 eV of the Fermi energy. Its high charge 
carrier mobility and peculiar electronic properties have 
stimulated considerable research into the possibilities of 
using graphene for electronic and spintronic applications. 

Graphene is often treated theoretically as a free- 
standing two-dimensional sheet. Though this often ap- 
pears to be a reasonable model for describing observed 
properties, in many experimental situations there is some 
deviation from this ideal because there is some form of 



physical contact with the environment. This can consist 
of atomic and molecular impurities in or on the graphene 
sheet, contact with an ins ulating sub strate, a gate elec- 
trode or metallic leads, etc. 12121151161121 ■\Yhile the Fermi en- 
ergy of free-standing graphene coincides with the conical 
points, adsorption on subst rates can alter its ele ctronic 
properties significantly. i9, , 20|2i|22 | 23|24|25|26 | 27|28| 

ample, the weak interaction of graphene adsorbed on the 
(0001) surface of insulating hexagonal boron nitride (h- 
BN) is enough to destroy graphene's characteristic con- 
ical points and open a band gap of some 50 meV^ 
Even when the interaction is sufficiently weak to leave 
the conical points essentially unchanged, it can lead to 
a larg e shif t of the Fermi energy away from the conical 
points.E^ll] 

Since measurement of the electronic transport prop- 
erties ofgra^henere^uires making contacts with metal 
leadsl ^'^'^i^^ ' ^^'^"''^^ ' ^^ ' '^^ ! it is important to understand 
such electronic and structural properties as the charge 
transfer between graphene and the metal substrate, the 
graphene-metal binding energies, distances etc. Charge 
transfer at a metal-graphene interface results in doping of 
the graphene sheet. Because the sign and the magnitude 
of the doping depend upon the metal, p-n junctions can 
be realiz ed by attaching electro des of different metals to 
graphcnc .'^'^'^^"^^'^^'^'^"^^'^^"^"'^^' 

There have been numerous theoretical and experimen- 
tal studies on semiconducting carbon nanotubes con- 
ta cted to m etals such as Al, Au, Pt, Pd, Ca and 
rji^ |42|43 | 4 4i45 46 gjj-^gg ^ graphene sheet can be considered 
as a carbon nanotube of infinite radius, the chemical in- 
teraction between graphene and metal substrates can be 
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FIG. 1: (Color online) (a) The most stable symmetric configuration of graphene on Cu, Ni and Co (111) lias one carbon atom 
on top of a metal atom (A site), and the second carbon on a hollow site (C site), (b) Graphene on Al, An, Pd and Pt(lll) 
can be modeled in a 2 x 2 graphene supercell with 8 carbon atoms and 3 metal atoms per layer. Shown is the most stable 
symmetric geometry. 



expected to be similar to that between metal contacts 
and nanotubes. 

In this paper we use first-principles calculations at the 
level of density functional theory (DFT) to characterize 
the adsorption of graphene on a variety of metal sub- 
strates. A preliminary account of our results was given 
in Ref. Hi The (111) surfaces of Al, Co, Ni, Cu, Pd, 
Ag, Pt, Au and the Ti(OOOl) surface cover a wide range 
of work functions and different types of chemical bond- 
ing, which allows for a systematic study of the metal- 
graphene interface. We focus on the interaction and 
charge transfer between graphene and the metal sub- 
strate, and in particular the effects they have on the 
doping of graphene by the metal. Because the charge re- 
distribution at the graphene-metal interface can be char- 
acterized experimentally by measuring the work function 
of the graphene-covered metal, we also calculate the work 
functions of these systems. 

The structural details of the metal-graphene interfaces 
will be presented elsewhere.'^ The most important re- 
sult for the purposes of this study is that there are two 
classes of graphene-metal interfaces. Whereas graphene 
is chemisorbed on Co, Ni, Pd and Ti, the binding to 
Al, Cu, Ag, Au and Pt is much weaker. The electronic 
structure of graphene is strongly perturbed by chemisorp- 
tion but is essentially preserved in the weak binding "ph- 
ysisorption" regime. For physisorbed graphene there is 
generally electron transfer to (from) the metal substrate, 
causing the Fermi level to move downward (upward) from 
the graphene conical points. This can be viewed as dop- 
ing graphene with holes (electrons) by adsorption. 

Naively one might expect the type and amount of 
doping to depend only on the difference between the 
work functions of free-standing graphene and of the clean 
metal surface. At typical equilibrium separations, the 



potential profile and therefore the doping are, however, 
altered significantly by an interface dipole arising from a 
direct short-range metal-graphene interaction. Using the 
DFT results, we develop an analytical model that quan- 
titatively describes the doping of physisorbed graphene. 
This model also predicts how physisorption of graphene 
modifies the metal work function. 

In order to characterize the doping of chemisorbed 
graphene a different approach must be used. Since 
chemisorption perturbs the electronic structure of 
graphene strongly, doping cannot be simply deduced 
from the shift of the Fermi level with respect to the 
conical points. Instead, we consider the work function 
of the graphene-covered metal, which is always a well- 
defined quantity. In a current-in-plane transport exper- 
iment only part of the graphene sheet covers (or is cov- 
ered by) the metal electrode, whereas an adjacent part 
is free-standing. The difference between the work func- 
tion of the graphene-covered metal electrode and free- 
standing graphene then determines the direction of the 
charge transfer between these two parts and hence the 
doping. According to this model, graphene is doped n- 
type by Co, Ni, Pd and Ti contacts. 

The paper is organized as follows. In Sec.|lT]we state 
the most important computational details of the den- 
sity functional calculations and summarize in Sec. [Ill A| 
the key results of a more extensive study of the bi nding 
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of graphene to various metal substrates!^ Section 
contains results of the first-principles calculations for 
the doping and work function of graphene adsorbed on 
these different substrates. A phenomenological model to 
describe the doping and work function of physisorbed 
graphene is introduced in Sec. |III C| and in Sec. |IIID| 
chemisorbed graphene is discussed. The sensitivity of 
the results to the computational approximations used is 
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TABLE I: and a°^^ represent the experimental cell parameters of the surface unit cells shown in Fig. ^s^) and Fig. [ijb) , 
respectively, for graphene on various metals. All calculations are performed with the lattice constant of graphene optimized 
using the LDA, ahox = 2.445 A. The calculated equilibrium separation doq is the separation in the z direction between the 
carbon atoms of the graphene sheet and the relaxed positions of the topmost metal layer, averaged where applicable over the 
carbon and metal atoms in the lateral supercell. The binding energy AE is the energy per carbon atom required to remove the 
graphene sheet from the metal surface. Wm and W are, respectively, the calculated work functions of the clean metal surfaces 
and of free-standing and adsorbed graphene, and W^'^, W°^^ are the corresponding experimental values. AEf is the Fermi 
level shift of physisorbed graphene. 
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discussed in Sec. lIIIEl A short discussion and conclusions 



are presented in Sec. IV 



II. COMPUTATIONAL DETAILS 

We calculate DFT ground state energies and optimized 
geometries using a plane wave basis set and the PAW 
formalism at the level of the local (spin) density ap- 
proximation, L(S)DA,^^ as implemented in the VASP 
code.^'^ ^'"^ '''^ The plane wave kinetic energy cutoff is 
set at 400 eV. A metal surface is modelled in a supercell 
as a finite number of layers of metal plus a region of vac- 
uum repeated periodically in the direction perpendicular 
to the layers. The supercell used to model the graphene 
metal adsorption is constructed from a slab of six layers 
of metal atoms with a graphene sheet adsorbed on one 
side and a vacuum region of ~ 12 A. A dipole correction 
is applied to avoid spurious interactions between periodic 
images of the slab.^^ 

We choose the in-plane lattice constant of graphene 
equal to its optimized LDA value, a = 2.445 A, adapt- 
ing the lattice constants of the metals accordingly. The 
graphene honeycomb lattice then matches the triangular 
lattice of the metal (111) surfaces in the lateral unit cells 
shown in Fig.[T| The approximation made by this match- 
ing procedure is reasonable since the mismatch with the 
lattice parameters of the metal (111) surfaces is only 0.8- 
3.8%, as seen in Table [ij In optimizing the geometry, the 
positions of the carbon atoms as well as those of the top 
two layers of metal atoms are allowed to relax. All re- 
sults reported in this paper are obtained for structures 
adapted to the LDA optimized in-plane lattice constant 



of graphene. 

We use the tetrahedron schem^^for accurate Brillouin 
Zone (BZ) integrations, sampling the BZ of the small 
and large cells in Fig. [T] with 36 x 36 and 24 x 24 k- 
point grids, respectively, and explicitly including the F, 
K and M high symmetry points. Note that on doubling 
the graphene lattice vectors to match those of Au, Pt, 
Cu, Ag, Al and Pd, the K point corresponding to the 
primitive unit cell of graphene is folded down onto the 
K point of the smaller Brillouin zone. The electronic 
self-consistency criterion is set to 10~^eV. Such a strict 
convergence is required to obtain accurate forces, which 
are essential in order to obtain reliable optimized struc- 
tures. Total energies are converged to within 10^^ eV in 
respect of ionic relaxation. Explicit total energy calcu- 
lations show that the structures in Fig. [T] represent the 
most stable symmetric configurations of graphene on the 
metal substrates studied, in agreement with experimen- 
tal results where available .li^ 

Detailed interfaces structures will be reported 
elsewhere.!^ Here we note that the L(S)DA functional 
gives a much better description of graphene-metal sub- 
strate binding energies and equilibrium distances than 
the commonly used generalized gradient approximation 
(GGA) functionals. Since the work functions, calculated 
with the L(S)DA, of clean metal surfaces and of those 
covered with graphene are sufficiently accurate, we use 
the L(S)DA functional. In Sec. IV we will show that, 
provided the graphene-metal substrate equilibrium sepa- 
ration is obtained correctly, the charge transfer and con- 
sequently the doping of graphene do not depend strongly 
on the choice of density functional. 

Matching the graphene lattice with the Ti(OOOl) sur- 
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FIG. 2: (Color online) Band structures of graphene adsorbed upon An, Pt, Cu, Ag, Al, Pd, Ni and Co (111) substrates. The 
Fermi level is at zero energy. The amount of carbon pz character is indicated by the blackness of the bands. The conical point 
corresponds to the crossing of bands at K with predominantly pz character, as is clearly visible for (physisorbed) graphene on 
Au, Pt, Cu, Ag, and Al. For (chemisorbed) graphene on Pd, Ni, and Co, the conical points disappear and the bands have a 
mixed character. The labels MIN/MAJ indicate the majority and minority spin bands of graphene on Ni and Co. The first 
and third top panels correspond to the band structure of free-standing graphene calculated with the primitive unit cell and 
2x2 graphene supercell, respectively. Inset: the two-dimensional Brillouin zones of graphene for the primitive unit cell and 
2x2 graphene supercell. In the supercell the bands are downfolded, the area enclosed by the bold lines in the primitive BZ 
translates to the corresponding one in the supercell BZ. 



face is more difficult since there is a lattice mismatch of 
20%. To accommodate this mismatch we use a graphene 
7x7 lateral supercell and a BZ sampling of a similar 
density as above. The equilibrium separation given in 
Table |T]is the value obtained by averaging over the lateral 
supercell. The details of the graphene/Ti(0001) calcula- 
tions will be reported elsewhere. Here we focus on the 
charge redistribution at the interface and the doping of 
graphene. 



III. RESULTS 

A. Metal-graphene binding 

The calculated equilibrium bonding distances, the 
binding energies and the work functions for adsorption 
of graphene on all metal substrates studied in this paper 
are listed in Table |l] The binding energies Ai? and equi- 
librium separations doq immediately show that the metals 
can be divided into two classes. For graphene adsorbed 
on Co, Ni, Pd(lll) and Ti(OOOl), I^E > 0.1 eV/carbon 



atom and dgq ^ 2.3 A. In contrast, adsorption on Al, 
Cu, Ag, Au and Pt(lll) leads to much weaker bonding, 
AE < 0.04 eV/carbon atom, and larger equilibrium sep- 
arations, (icq ^ 3.3 A. The equilibrium geometries and 
distances obtained are in agree ment with available ex- 
perimental data and calculationJ^SHnHHMI and appear to 
be similar to the bonding found between graphene and 
carbon nanotubes; carbon nanotubes are usually bonded 
strongly to Pd and Ti w hereas the bo nding with Al, Ag, 
Au, Ca and Pt is weaker .I^SSIIIIIIIIII 

The difference between the two classes of metal sub- 
strates is reflected in the electronic structure of adsorbed 
graphene as shown in Fig. [2j When the binding energy 
is large, i.e., if graphene is adsorbed on Co, Ni, Pd or 
Ti, the graphene bands are strongly perturbed. In par- 
ticular, the characteristic conical points of graphene at 
K are destroyed. Graphene Pz-states hybridize strongly 
with the metal d-states and the corresponding bands ac- 
quire a mixed graphene-metal character. It demonstrates 
that graphene is chemisorbed on these substrates. 

In contrast, if the metal-graphene interaction is 
weaker, i.e., when graphene is adsorbed on Al, Cu, Ag, 
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FIG. 3: (Color online) Schematic illustration of the parame- 
ters used in modelling the interface dipole and potential step 
formation at the graphene-metal interface. 



Au or Pt, the graphene bands, including their conical 
points at K, can still be clearly identified. We reserve 
the term physisorption to describe this type of bonding. 
Unlike in the case of free-standing graphene where the 
Fermi level coincides with the conical point, physisorp- 
tion generally shifts the Fermi level. Even when there is 
no interaction or the interaction is weak, this does not 
preclude the transfer or charge between graphene and 
the metal substrate resulting from the equilibration of 
the chemical potentials. 



B. Doping of physisorbed graphene 

In physisorbed graphene the conical points in the 
graphene band structure are preserved, but charge trans- 
fer to or from the metal substrate shifts the Fermi level. 
A schematic representation of the parameters we use to 
describe this situation is shown in Fig. |3]for the case of 
electron transfer from graphene to the metal. A shift 
upwards (downwards) with respect to the conical points 
means that electrons (holes) are donated by the metal 
to graphene, making the latter n-type (p-type) doped. 
We extract the Fermi level shifts AEp of graphene ph- 
ysisorbed on a number of metals from the band structures 
shown in Fig. [2] and plot them in Fig. [4] At equilibrium 
separations from the metal substrates, graphene is doped 
n-type on Al, Ag and Cu, and p-type on Au and Pt. 
In the following section we develop a phenomenological 
model to describe these first-principles results. In the 
remaining part of this section we identify the physical 
parameters that play a role in this model. 

The work function W of a, graphene-covered metal is 
given by the position of the Fermi level {W = —Ef). Be- 
cause the density of states of graphene is so small com- 
pared to that of the local density of states of a typical 
transition metal surface, the shifts required to equilibrate 



the Fermi levels when charge transfer occurs take place 
almost entirely in graphene. For physisorbed graphene 
where the interaction is so weak that its electronic struc- 
ture is unchanged, W should be related to the Fermi level 
shift in a simple way 



AEf = W-Wg, 



(1) 



where Wq is the work function of free-standing graphene. 
The work function shifts are calculated separately and 
are plotted in Fig. |4] We see that while Eq. ([TJ holds 
for a relatively large separation d — 5.0 A between the 
graphene sheet and the metal surface, there is a small 
deviation of ^ 0.08 eV at the equilibrium separation 
d « 3.3 A which can be traced to a perturbation of the 
graphene electronic structure by physisorption, that can- 
not be described as a rigid shift. In the following discus- 
sions we will neglect this small (non-rigid shift) pertur- 
bation. When graphene is chemisorbed and the Fermi 
level shift cannot be determined from the strongly per- 
turbed band structure, the work function W is still a 
well-defined parameter. 

Because the work functions of graphene, Wq, and of 
most metal surfaces, Wmi differ, electrons are transferred 
from one to the other to equilibrate the Fermi levels if 
the two systems communicate. Charge transfer between 
metal and graphene results in the formation of an inter- 
face dipole layer and its associated potential step, AV. 
We can use the plane-averaged electron densities n{z) to 
visualize the electron redistribution upon formation of 
the interface 



An{z) = nM\G{z) - nuiz) - naiz), 



(2) 



where nMioiz), nuiz) and nQ{z) denote the plane aver- 
aged densities of the graphene-covered metal, the clean 
metal surface, and free-standing graphene, respectively. 
Notice that the structure of the clean metal surface is 
required to be the same as that of the graphene-covered 
metal surface. The results for graphene physisorbed on 
Al, Cu, Ag, Au and Pt, are shown in Fig. [5] An is lo- 
calized near the interface for all metal substrates and in 
the majority of cases it has the shape of a simple dipolar 
charge distribution. 

We estimate the charge q (per carbon atom) that is 
responsible for the dipole by integrating An from the 
node at zq between the metal surface and the graphene 
sheet, 



dz An(z) /Nq 



(3) 



where iVc is the number of carbon atoms in the unit cell; 
— e is the charge of an electron. These numbers are in- 
cluded in Fig. |5] The sign and size of the dipole charges 
are consistent with the changes of the metal work func- 
tion upon adsorption of graphene. Note that relatively 
small values of charge transfer give rise to quite substan- 
tial work function changes, see Table [ij 
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FIG. 4: (Color online) Calculated Fermi energy shift with re- 
spect to the conical point, Ai^p (dots), and W — Wa (crosses) 
as a function of the clean metal-graphene work function differ- 
ence Wm — VKg- The lower (black) and the upper (green/grey) 
points are for the equilibrium (~ 3.3 A) and large (5.0 A) 
graphene-metal surface distances, respectively. The solid and 
dashed lines follow from the model of Eq. ([7|. The insets 
illustrate the position of the Fermi level with respect to the 
conical point. 



The above analysis points to the use of a plane capaci- 
tor model to describe the potential step AV. As sketched 
in Fig. [Sj the charge distribution is then modelled as two 
sheets of charge ±q. Since the charge is predominantly 
localized between graphene and the metal surface, the 
effective distance Zd between the charge sheets should be 
smaller than the graphene metal separation d. 

If the interaction between graphene and the metal sur- 
face is weak, as in the case of physisorption, one naively 
expects that electrons will be transferred to graphene 
if the clean metal work function is lower than that of 
free-standing graphene, i.e., if Wm < Wq- Electrons 
should then flow from graphene to the metal surface if 
Wm > Wq and the crossover point from n- to p-type dop- 
ing would be exactly at Wm = Wg ■ The results obtained 
for the equilibrium graphene-metal separation, d ~ 3.3 A 
in Fig. [4] do not confirm this simple picture of the doping 
mechanism. Instead, the crossover point between n- and 
p-type doping is found for a metal with a work function 
Wm = Wg + 0.9 eV. 

This simple picture of charge transfer cannot be en- 
tirely wrong. If the graphene-metal separation is in- 
creased, the crossover point from n- to p-type doping 
decreases to its expected value, Wm ~ Wq, for large sep- 
arations. This is illustrated by the upper curve in Fig. l4] 
calculated for a graphene-metal separation of d — 5.0 A7 
It clearly indicates that, at the equilibrium separation 
deq ~ 3.3 A, the charge reordering at the graphene-metal 
interface is the result not only of a charge transfer be- 
tween metal and graphene electronic levels that equili- 
brates the graphene and metal Fermi energies but that 
there is also a contribution from a direct interaction be- 
tween the metal and graphene. A similar interaction. 
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FIG. 5: (Color online) Plane-averaged electron difference den- 
sity An{z) (per unit cell) showing the charge displacement 
upon physisorption of graphene on M(lll) surfaces where M 
— AI, Ag, Cu, Au, and Pt. g/e is the number of electrons per 
carbon atom calculated by integrating An{z) from the central 
node to infinity.'^ 



which has a significant repulsive contribution, plays an 
important role in describing the dipole formation for 
closed shell at oms a nd organic molecules adsorbed upon 
metal surfaces .^^^l^ The interaction depends on the wave 
function overlap between the metal and the adsorbed 
species. We expect it therefore to be very sensitive to 
the metal-graphene separation d and to vanish exponen- 
tially with increasing d. 



C. Phenomenological model 

In this section we construct a simple and general model 
to describe the Fermi level and work function shifts cal- 
culated from first-principles for graphene physisorbed 
on AI, Ag, Cu, Au, and Pt substrates. All relevant 
parameters are shown in Fig. [3] We start by writ- 
ing the work function of the graphene-covered metal as 
W{d) = Wm — A.V{d), where AV is the potential step 
generated by the interface dipole layer. Its size depends 
on the graphene-metal separation d. The Fermi level 
shift in graphene and the work function are related by 
Eq. ([ij. Use of these relations implicitly assumes that 
the graphene electronic energy levels around the Fermi 
energy are essentially unchanged by the interaction be- 
tween graphene and the metal and that the band struc- 
ture of graphene is just rigidly shifted by the interface 
potential AV. 

A key element of the model is to write the interface 
potential step as AV{d) = Atj:{d)+Ac{d). The first term, 
Atr(d), results from the direct charge transfer between 
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graphene and the metal, which is driven by the difference 
in work functions. The second term, Ac(d), describes the 
short-range interaction discussed in the previous section, 
that results from the overlap of the metal and graphene 
wave functions. We parameterize it as 



Ac(d) = e"^"^ (ao + aid + ajd^). 



(4) 



i.e., we assume that it vanishes exponentially with in- 
creasing graphene-metal separation d. The exact asymp- 
totic functional dependence of Ac(d) for large d would 
very likely change if one were to go beyond DFT/LDA 
and take the van der Waals interaction into account. The 
asymptotic form of Ac(d) is, however, not important be- 
cause |Ac(d)| becomes negligible for large d anyway. 

To model the electron transfer contribution, Atr((i), 
we use a plane capacitor model so Atr(rf) = aN{d)zd, 
where a = e^je^A = 34.93 eV/A with A = 5.18 
the area of the graphene unit cell, and is the num- 

ber of electrons (per unit cell) transferred from graphene 
to the metal. Note that A^(fi) becomes negative if elec- 
trons are transferred from the metal to graphene. The 
parameter is the effective distance between the sheets 
of transferred charge on graphene and the metal. It is 
smaller than the geometrical metal-graphene separation, 
Zd < d, because most of the charge is localized between 
the graphene layer and the metal surface. We approxi- 
mate the effective distance between the charge sheets by 
Zd = d — do with do a constant. 

A closed set of equations for the Fermi level shift AE-p 
and the work function W is obtained by determining the 
relation between AEp and the number of electrons N 
transferred between graphene and the metal. For an 
energy range within ±1 eV of the conical points, the 
graphene density of states is described well by a linear 
function 



D{E) = Do\E\, 



(5) 



with Dq = 0.09/(eV^ unit cell). Integrating the density 
of states from the neutrality point, dE D{E), yields 

the required relation, N = sign{AEY)DQAEl /2P^ 
The model can be summarized by the set of equations 



'W{d) = WM-AV{d), 
AF(d) = At,(d) + A,(d), 
AtM)^aN{d) (d-do), 
N{d) = sign(Ai;F) ^DoAEpidy 

^E^id) = W{d) - Wq. 



(6) 



Solving this set results in the following simple expression 
for the Fermi level shift 



^Vl + 2aDo{d-do)\WM-WG-A,{d}\-'^ 

^^^('^) = ^ aDo{d-do) ' 

(7) 

where the sign of AE-p is determined by the sign of Wm — 
Wq — A(.. The work function of the graphene-covered 
metal surface is then obtained from Eq. Q. 
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FIG. 6: (Color online) Top panel: Fermi level shifts relative to 
the Dirac point, Aii'F(d), as a function of the graphene-metal 
separation d for physisorbed graphene. Bottom panel: calcu- 
lated work functions W{d) relative to that of a free graphene 
sheet, Wg- The dots (top) and crosses (bottom) give the cal- 
culated DFT/LDA results, the solid lines describe the results 
obtained from the phenomenological model, Eq. ItI.'^ 



We used first-principles calculations to obtain AEp 
and W explicitly for a range of separations d. These first- 
principles data points are shown in Fig.[6j We could now 
fit the model to these points and so obtain a simple inter- 
pretation of the numerical results. However, we can do 
better than that. It turns out that the parameter do and 
the function Ac{d) depend only very weakly on the metal 
substrate. This means that we can fit these quantities to 
the first-principles results for a single metal substrate and 
subsequently use Ac{d) and do as universal parameters 
to predict the Fermi level shifts in graphene for all metal 
substrates. To determine do and Ac(d) we use the DFT 
results for graphene on the Cu(lll) surface.'^ 

The Fermi level shift for graphene physisorbed on any 
other metal substrate can then be obtained from the 
model, using only the work function of the clean metal 
surface, Wm, and that of free-standing graphene, Wq, as 
input parameters in Eq. Q . The accuracy of this model 
is demonstrated by Fig. ^and the top panel of Fig. |6] 
The latter shows that the distance dependence of the 
Fermi level shift is represented very well by the model. 
For graphene on Cu, Ag, and Au(lll) the deviations 
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from the first-principles results are <C 0.1 eV, whereas for 
graphene on Al and Pt(lll) they are still < 0.1 eV. The 
slightly larger deviations for Al and Pt might be due to 
a more complex interaction between graphene and these 
surfaces, caused by the open p and d shells of Al and Pt, 
respectively. The latter might also be responsible for a 
non-monotonic behavior of A£'F(d) and W{d) for Al and 
Pt in Fig. [§ 

Once Ai?F has been determined, the work function 
W of the metal-graphene system can be obtained using 
Eq. pj and the results are shown in the bottom panel of 
Fig. [or The work function given by the model agrees with 
the first-principles results within < 0.2 eV for d > 3.3 A, 
the equilibrium separation, but seems to be systemati- 
cally lower for smaller values of d. The difference be- 
comes smaller upon increasing d. For small graphene 
metal separations d, Eq. ([T]) no longer holds exactly, as 
the electronic structure of graphene is perturbed by the 
presence of the metal substrate. Note, however, that 0.2 
eV is comparable to the difference between the experi- 
mental and the DFT/LDA work functions, see Table |l) 
which means that the error resulting from the model is 
tolerable. 

According to Eq. (|6| one can calculate the sign and 
concentration of the charge carriers in graphene, N, 
from the Fermi level shift Ai?p. The nominal charge on 
graphene per unit cell containing two carbon atoms, eA^, 
is —14.6, —1.6, —1.3, 1.6, 4.9 x lO^'^e for graphene on 
Al, Ag, Cu, Au, Pt(lll), respectively.'^! Although these 
charges are very small, the Fermi level shifts are quite 
substantial because the graphene DOS close to the con- 
ical points is so low, but they are still within the linear 
regime described by Eq. (|5]). In the terminology used in 
semiconductor physics this amount of charge per unit cell 
would be classified as heavy doping. 

The crossover from p- to n-type doping occurs when 
the Fermi level coincides exactly with the conical points 
of graphene, i.e., AE'p = 0. According to Eq. ([t]), this 
happens if the work function of the metal is given by the 
critical value 



Wo(d) = W^G + Ac(d). 



(8) 



The critical work function Wq depends on the graphene- 
metal separation d, since the term Ac resulting from 
the direct graphene-metal interaction depends strongly 
on d, becoming negligible if d ^ 4.2 A. Wq then ap- 
proaches Wq = 4.5 eV, i.e., the critical work function is 
that of free-standing graphene. However, at the equihb- 
rium graphene metal separation, d^q = 3.3 A, A^, 0.9 
eV, leading to a critical work function WQ{d) « 5.4 eV. 
Though the short-range graphene-metal interaction does 
not significantly change the graphene band structure, it 
does lead to a sizeable potential step at the equilibrium 
separation which is downwards from metal to graphene as 
indicated in Fig. [3] The size of this potential step is rela- 
tively insensitive to the metal substrate. A similar poten- 
tial step has been observed in the adsorption of closed- 
shell molecules on metal surfaces, where it has been in- 





7.5 




5.0 




2.5 







o 


-2.5 


X 


-5.0 




-7.5 




-1.0 







q=0.013e 





5 10 15 20 

z(A) 





- 


q=0.021e| 


> 



- Co J 




_ q=0.022e 


> 



10 15 

z(A) 



20 



10 15 

z(A) 



20 



FIG. 7; (Color online) Plane-averaged electron difference den- 
sity An(z) (per unit cell) showing the charge displacement 
upon formation of the chemisorbed graphene-M(lll) inter- 
face, with M = Ni, Co, Pd. q is the charge (per carbon atom) 
calculated by integrating An(2) from the central nodal point 
to infinity.'^ 



terpreted in terms of an exchange repulsion between the 
electrons on the molecules and the metal substrate.!^ 

The phenomenological model we have outlined de- 
scribes the doping of graphene by metal contacts and the 
work function shifts caused by adsorption of graphene 
when the graphene-metal bonding is weak. The model 
is based on the linearity of the graphene DOS, which 
holds for energies within ±1 eV around the conical points. 
Therefore, the criterion for the validity of the model is 
Ag ^ |Ai?F| ^ 1 eV, where Ag is a band gap induced in 
graphene by interaction with the substrate.!^ 



D. Chemisorbed graphene 

In the previous sections we defined the doping of ph- 
ysisorbed graphene in terms of the Fermi level shift AiJp 
with respect to the conical points in the graphene band 
structure. Negative and positive AiJp correspond, re- 
spectively, to n-type and p-type doping. This procedure 
cannot be used for graphene that is chemisorbed on the 
Ni, Co, or Pd(lll) or on the Ti(OOOl) surface because 
the strong graphene-metal bonding interaction destroys 
the conical points; see the Ni, Co and Pd panels in Fig.|2] 

The more complex bonding of chemisorbed graphene is 
also illustrated by comparing the plane-averaged electron 
difference densities An(z) shown in Figs. [5] and [7j In the 
physisorption case. An still has the characteristics of a 
simple interface dipole while in the chemisorption case 
it is much more complicated, indicating the formation 
of new bonds at the interface between graphene and Ni, 
Co and Pd. The charge reordering at the interface upon 
chemisorption is substantial, as reflected by the larger 
values for q given in Fig. |7] This leads to considerable 
shifts in the metal work functions upon chemisorption of 
graphene. In all cases studied in this paper graphene acts 
as an electron donor, lowering the metal work function. 
For Ni, Co and Pd(lll) the work function lowering is 
1.81, 1.66, and 1.64 eV, respectively, and for Ti(OOOl) it 
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is 0.56 eV. 

Despite the impossibility of identifying a Fermi level 
shift and therefore the type and magnitude of the dop- 
ing from a simple examination of the band structures 
in chemisorbed graphene, it is possible to define a suit- 
able measure for application to current-in-plane (CIP) 
geometries. In a CIP geometry only part of the graphene 
sheet is covered by metal electrodes, and the greater part 
of the sheet is "free-standing". The type and effective 
concentration of charge carriers in graphene contacted 
to metallic leads can be measured in experiments us- 
ing the CIP geometry shown schematically in Fig. [s] (for 
clarity, for the physisorbed case). At a large distance 
from the metal contact the Fermi level in free-standing 
graphene approaches the conical points. At the metal- 
lic contact the Fermi level is fixed by the interaction 
with the metal electrode. The difference between the 
Fermi levels in the adsorbed and free-standing graphene 
is given by the difference between the work function of 
the graphene-covered metal surface, W, and that of free- 
standing graphene, Wq, or in other words, by Eq. Q; see 
Fig. ^a). We already used this relation in our descrip- 
tion of physisorbed graphene. Since the work function W 
can be determined for chemisorbed as well as physisorbed 
graphene, it can be applied to all metal electrodes in the 
CIP geometry. 

To accommodate the Fermi level difference, charge 
transfer takes place between the contacted and free- 
standing regions.*^ Alignment of the Fermi levels results 
in the band bending shown schematically in Fig. [sjb). 
The amount of band bending is given by Eq. ^. The 
band bending region is p- or n-type doped, depending 
on the sign of AE'p. If AEp > 0, graphene is p-type 
doped, and if AEp < 0, it is n-type doped. Different 
metal electrodes can then be used to make p-n junctions 
in graphene. 

For graphene chemisorbed on Ni, Co, Pd(lll) and 
Ti(OOOl) surfaces, we find AEp = W - Wq ^ -0.82, 
-0.70, -0.45, and -0.34 eV, respectively, see Table |l] 
Chemisorption of graphene on these surfaces lowers their 
work function to below that of free standing graphene. 
Therefore, we expect graphene to be n-doped by these 
metal electrodes. 



E. Sensitivity to approximations 

The calculations on the interaction between graphene 
and the metal substrates we have discussed so far, are at 
the level of DFT/L(S)DA. The results given in Table |l] 
show that calculated work functions of clean metal sur- 
faces agree with experimental data within ^ 0.2 eV. A 
similar agreement is observed between calculated work 
functions of graphene-covered metal surfaces and avail- 
able experimental data. The change in work function 
upon graphene adsorption is determined by the forma- 
tion of an interface dipole and the charge transfer be- 
tween metal and graphene. Apparently this charge re- 



vacuum level 




FIG. 8: (Color online) Schematic representation of a graphene 
sheet partly in contact with a metal electrode. The work 
function of the graphene-covered metal electrode, W , is here 
smaller than the work function of graphene, Wg, and the 
graphene sheet becomes n-type doped. Far from the elec- 
trode the conical point of graphene (bold dots) approaches 
the Fermi energy, (a) In the absence of communication be- 
tween the graphene covered metal and the intrinsic graphene 
sheet, there is a discontinuity AEy of the Fermi energies, (b) 
As soon as the two systems can communicate, equilibration 
of the Fermi energies takes place by the transfer of electrons 
from the low to the high work function system and the joint 
Fermi energy is fixed by the graphene-covered metal electrode, 
Ef = —W. This rearrangement of charge gives rise to a po- 
tential shift AV' . The band-bending and graphene doping 
depend on the distance x from the contact. 



distribution at the interface is described rather well by 
the LDA functional. 

Semi-local GGA junctionals are frequently preferred 
in DFT calculations. For graphene that is physisorbed 
on a metal surface the commonly used GGA function- 
als give an interaction that is either too weak or even 
purely repulsive making it impossible to predict the equi- 
librium distance between graphene and the metal sur- 
face. In order to test whether the Fermi level shift in 
adsorbed graphene and its work function are sensitive to 
the particular choice of the density functional, we have 
calculated these parameters for graphene on Cu(lll) as 
a function of the gr aphene- met al surface distance, using 
the PW91 GGA functional.l^^ The results are shown 
in Fig. [9] The Fermi level shifts AEp calculated with 
GGA are within ^ 0.07 eV of the ones obtained with 
LDA. The same holds for the work function difference 
W — Wgi which means that the description of the dop- 
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FIG. 9: (Color online) The dashed lines correspond to the 
interpolation curves for AEf and W{d) — Wg as a function 
of the graphene-Cu(lll) surface distance d calculated with 
the PW91 GGA functional. The dash-dotted lines represent 
the results obtained with the LDA functional with graphene 
stretched to the Cu(lll) lattice constant ahox = 2.49 A. As a 
reference, the solid lines give the LDA results obtained with 
the optimized graphene lattice constant Ohcx = 2.445 A. 



ing of graphene is not sensitive to the functional used in 
the calculations. However, the absolute values of W and 
Wg calculated using the GGA differ from those obtained 
using the LDA by ~ 0.1 — 0.3 eV, which is a typical dif- 
ference between work functions obtained with these two 
density functionals.!^ 

Another issue is the possible effect of a truly non-local 
van der Waals interaction between the graphene sheet 
and metal surface, which is neglected in our study. Van 
der Waals interactions have been addressed recently in 
calculations on graphite, hexagonal boron-nitride and di- 
atomic molecules of inert gases.^SMi A non-local correc- 
tion to a GGA-type density functional has been proposed 
in Ref. 69 that correctly reproduces the asymptotic van 
der Waals tail of the binding energy curve at large inter- 
molecular distances. This non-local correction has little 
impact on the charge distribution, however.'^ Since it is 
the charge distribution that gives the Fermi level shift, 
the work function and the doping of adsorbed graphene, 
we conclude that these quantities are adequately de- 
scribed by local or semi-local functionals. 

In our calculations we choose the in-plane lattice con- 
stant of graphene equal to its optimized LDA value 



a = 2.445 A, and adapt the lattice constants of the met- 
als accordingly. The approximation made by this match- 
ing procedure seems reasonable, since the mismatch with 
the lattice parameters of the metal (111) surfaces is only 
0.8-3.8%, see Table [ij The largest lattice mismatch is 
that between graphene and Cu(lll). To estimate the 
error on the Fermi level shift and the work function of 
adsorbed graphene, we have calculated these quantities 
while stretching graphene to the LDA optimized in-plane 
lattice constant of Cu(lll) (2.49 A). The results are 
shown in Fig. [9] The Fermi level shifts and work function 
differences are within 0.15 eV of the results obtained with 
the optimized graphene lattice constant. One expects to 
see a smaller effect for other metals, as the lattice mis- 
match between graphene and other metals is smaller. 



IV. DISCUSSION AND CONCLUSIONS 

The theoretical study performed in the previous sec- 
tions assumes that the graphene sheet is adsorbed on a 
clean crystalline metal contact. We have used the model 
represented by Eq. ([t]) to describe the Fermi level shift 
in physisorbed graphene. Interpreting experiments that 
are not carried out in ultra-high vacuum requires some 
modifications because of impurities that will be present 
at the metal-graphene interface .1^^ The work function of 
the clean metal surface must then be replaced with that 
of the metal surface contaminated by water molecules, 
oxygen and/or nitrogen, for instance. The short-range 
graphene-metal surface interaction, represented by the 
potential term Ac, and the equilibrium separation dcq 
will also be affected by such adsorbates and the effects 
will depend critically upon the adsorbate concentration. 
The same obviously holds for chemisorbed graphene. A 
large concentration of adsorbates on the metal surface 
could break the chemical interaction between graphene 
and the metal. 

In some experiments the interface between the 
graphene sheet and metal surface might contain a thicker 
buffer layer consisting, for example, of water, or a metal- 
oxide.'^ That would certainly modify the graphene- 
substrate interaction. An input parameter to the model 
of Eq. ^ then is the work function of the metal with 
the buffer layer on top. Using the plane capacitor model 
one should replace a in Eq. (|7| by u/k, where k is the 
effective dielectric constant of the buffer layer. Obviously 
the distance Zd between the charge sheets has to be mod- 
ified accordingly to take the thickness of the buffer layer 
into account. In addition the potential term Ac, which 
represents the short-range interaction of graphene with 
the substrate should now reflect the graphene-buffer layer 
interaction. 

A metallic buffer layer consisting of, for instance, Ti 
is often used in experiments to establish a good contact 
between graphene and the electrodes. Such layers have a 
thickness of typically ^ 5 nm, which is sufficiently thick 
that the contact should be considered as a graphenc-Ti 
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contact. 

Using first-principles DFT calculations, we have sys- 
tematically studied the interaction and charge transfer 
between graphene and a range of metal surfaces with 
different work functions. We found that graphene is 
chemisorbed on Co, Ni, and Pd(lll) surfaces, and on 
the Ti(OOOl) surface and that this strong interaction per- 
turbs the electronic structure of graphene significantly. 
In contrast, adsorption of graphene on Al, Cu, Ag, Au 
and Pt(lll) surfaces leads to weak bonding, which pre- 
serves the typical graphene electronic structure, includ- 
ing the conical points. Even in this physisorbed case, 
however, there is a short-range graphene-metal interac- 
tion, as well as a charge transfer between the graphene 
and metal states. These result in a doping of graphene, 
i.e., a shift of the Fermi level with respect to the conical 
points, and a significant change in the work function of 
the graphene-covered metal surfaces, as compared to the 
clean metal surface. 

To extend the apphcability of the DFT results on ph- 
ysisorbed graphene we develop a simple general model, 
which takes into account the charge transfer between the 
graphene and metal states and the short-range graphene- 
metal interaction. We find that the latter only weakly 
depends on the metal. Therefore it can be fitted using 
the DFT results on one metal substrate. The model then 
only has the work functions of free-standing graphene and 
that of the clean metal surface as input parameters and it 
predicts the Fermi level shift and carrier concentration in 
graphene, as well as the work function of the graphene- 
covered metal substrate. We find that graphene is n-type 
doped if the metal work function Wm ^ 5.4 eV, whereas 



it is p-type doped if Wm ^ 5.4 eV. 

For the CIP geometry, where only part of the graphene 
sheet covers the metal electrode, we propose a definition 
of the doping that is based upon the work function of 
graphene-covered metal surface. This definition is valid 
both for the chemisorbed and physisorbed cases. It pre- 
dicts that adsorption of graphene on Al, Ag, Cu, Co, 
Ni, Pd(lll) surfaces and on the Ti(OOOl) surface leads 
to n-type doping. The high values of the work functions 
of Au and Pt(lll) substrates lead to p-type doping of 
graphene. 

Both the analytical model and the CIP definition of 
the graphene doping are derived from general principles 
and should be applicable to any metal surface on which 
graphene can be epitaxially grown. This opens up the 
possibility of a general understanding of p-n junctions 
prepared by doping graphene with metal contacts.^^ 
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